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Abstract 

It is shown that because of the radiation pressure a Schrodinger cat state 
can be generated in a resonator with oscillating wall. The optomechanical con- 
trol of quantum macroscopic coherence and its detection is taken into account 
introducing new cat states. The effects due to the environmental couplings 
with this nonlinear system are considered developing an operator perturbation 
procedure to solve the master equation for the field mode density operator. 



PACS number(s): 42.50.Dv, 42.50. Vk, 03.65.Bz 



On leave from Lebedev Physical Institute, Moscow 



1 



1 Introduction 



One of the fundamental aspects of quantum mechanics is the existence of inter- 
ference among quantum states which signs the difference between a superposition of 
states and a mixture of states. The quantum theory may adequately well describe 
macroscopic objects by means of a linear superposition of states with macroscopi- 
cally distinguishable properties. Recently, due to the improved technology, there was 
a growing interest on the possibility of observing such superposition states, commonly 
known as Schrodinger cats . Good candidates for these macroscopic states are the 
coherent states of an e.m. field mode. The properties of superposition of two generic 
coherent states has been studied in Ref. and the simplest superposition of even 
and odd coherent states was introduced in Ref. ^ . A review of these states is given 
in Ref. 0. Within the field of optics several proposals for the generation of linear 
superpositions of coherent states in various nonlinear processes || and in quan- 
tum non- demolition measurements [0 have been made. It is worthy to note that the 
field in a cat state has a lot of advantages in optical communication p. However, 
by coupling the system to its enviroment, as in the act of measurement, one always 
introduces dissipation and decoherence effects, which tend to destroy any quantum 
features |^. 

In the common scheme of the Kerr-like medium modeled by an anharmonic os- 
cillator, it was shown |TD| that the photon number distribution and interferences in 



phase space are highly sensitive to even small dissipative coupling. This fact, plus the 
smallness of the x^^^ nonlinearity, makes the prospect of experimentally producing 
and detecting such states highly questionable in these media. 

On the other hand, it is well known |ll]] that an empty optical cavity with a 
moving mirror may mimic a Kerr-like medium when it is illuminated with coherent 
light. The effect of intensity dependent optical path is due, in this case, to the 



radiation pressure force. 

In this paper we shall present such a model as an alternative one for the generation 
of Schrodinger cats. We will show that, with the appropriate measurement technique, 
it could also be useful for revealing quantum macroscopic coherence. 

2 The Model 

We consider a linear Fabry-Perot empty cavity with one fixed partially reflecting 
end mirror and one perfectly reflecting mirror, which can move (undergoing harmonic 
oscillations) under the influence of radiation pressure. If L is the equilibrium cavity 
length, the resonant frequency of the cavity will be 

Uc = TT^n (1) 

where n is an integer number determined by the frequency of the input light and c 
is the speed of light. We assume that the retardation effects, due to the oscillating 
mirror, in the intracavity field are negligible. We will also neglect the correction to 
the radiation pressure force due to the Doppler frequency shift of the photons [0. 
Thus we are able to write the Hamiltonian of the whole system as 

H = hujca^a + hiOrJi^h + Hint (2) 

where a, a'^'are the boson operators of the resonant cavity mode and 6, 6^ are the 
boson operators of the oscillating mirror with the mass m and the angular frequency 
ujm- This latter will be many order of magnitude smaller than uoc to ensure that the 
number of photons generated by the nonstationary Casimir effect [|T^ as consequence 



of the Casimir forces in the resonator with moving boundaries is completely 
negligible. Hint accounts for the fact that the intracavity photon changes its energy. 



by cUc, as the oscillating mirror moves ||15 



Hint = -hGa^a{h + 6^) (3) 



with the couphng constant given by 



/ * \ 1/2 



L \2'muJmj ^ ^ 

From the Hamiltonians (§) and @ we can derive, using the BCH formula for the Lie 



algebra []I6[, the time evolution operator in the following form 
where 

is the mirror quadrature operator, while 

E{t) = K^[t -sint]] F{t) = 2Ksm{t/2); K = G/um, (7) 

with t the time scaled by ujrn, i-e. we have replaced by t. Furthermore, from 
now on, we will consider the evolution operator omitting the free motion of the two 
modes a and b, i.e. the term inside the square brackets on Eq. (^. 

3 Generation of Schrodinger Cat States 

From Eq. (|^) one can immediately recognizes that the time evolution introduces 
anharmonicity due to the presence of the nonlinear term (a'^a)^ whose strenght de- 



pends also on time [0. It is also easy to see that at each time for which F{t) = 
the two subsystems are disentagled. Furthermore due to its macroscopicity we should 
consider the oscillating mirror initially in a thermal state at temperature T 

with z/{l — z) = Nth that represents the mean number of excitations of the mechanical 
oscillator, i.e. the number of thermal phonons. Thus starting from an initial coherent 



state |ao) for the radiation mode we have 



(9) 



with 



so that 



t* = 27rmi; mi e N 



F{t*) = 0; E{t*) = K^2nmi. 



(10) 



(11) 



Now in order to see the cat states, the fohowing condition must be fulfilled 



E{t*) = - + 27rm2; m2 E N 



(12) 



so that combining Eqs.(|Tl]) and (|T2|) one gets 



1 ,1 



K = — (- + m2) 



(13) 
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which can be read as a restriction on the possible values of the various external 
parameters. Thus if the above conditions are satisfied, we have 

1 r 



Pin = - [e-^'^/Vo) + e^^/^l - ao)\ [(-aole"^^/^ + {ao\e'''/^ 



(14) 



however, this is not the only way to create a quantum superposition in this system. 
In fact, let us consider the times t' for which 



TT 



E{t') = - + 271171; meN. 



(15) 



In these cases, obviously, F{t') is not necessarily zero then, the reconstruction of the 
superposed coherent states is impossible due to the entanglement between the two 
subsystems. One can now use a conditional measurement to create the desired states. 



performing a sort of quantum state engineering [T^ 
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Let us suppose that the mirror's quadrature x{t) is measured []T9|, giving the result 
Ut- The state of the radiation field after the measurement is found by projecting the 
system's state onto the eigenstate \yt) 

where C is a normalization constant 

C = i{yt\pT\yt))-' ■ (17) 
At the times t' we have, from Eqs. (|15|) and ( [T^ ) 

Pafterit') = ^ [c'^'^/Voe^^^*')^*' ) + C^^/^l - CoC^^^*')^*' ) 

which is a superposition of coherent states whose phase depends on the measurement 
process; and further, if the result of the measurement is 

y^' = lj^y (19) 



it is possible to recover in Eq. (^8|) the generalized even and odd coherent states 
like those discussed in [^, ^ which show quantum interference as well as other 
particular features. 

4 Quasiprobability and Marginal Distribution 

The evolved density operator of the whole system can be easily constructed by 
using the time evolution operator of eq. (^ 

pit) = Uit)\ao){ao\^PTU\t), (20) 

and then the evolution can be described for example, in terms of the Q-function 

oo 

Q(a,/5,t) = {a\{/3\pim\a) = e-H^-l"ol^-l'3P(l - z)Y.z^P\^ 

j=0 



n=0 



exp 



^E{t) - -F\t) 



ne 



it/2 



r!J(j-r)! 



(21) 



where the variables a, /3 are referred to the radiation and to the mirror respectively. 
However, since the distinguishing element of a linear superposition of coherent states 
is the presence of interference fringes in the marginal distribution, we are interested 
in that one, for the particular times discussed in the previous Section. Its definition, 
for a generic state pf^^^'^{t) of the radiation field, is given by 



P{X) = {X\pf^^'\t)\X) , 



(22) 



where \X) are eigenstates of the quadrature operator X = (a + a^)/2, while p^**^'*^ 
should be intended as Tr^jp} with Tr^ the trace over the mirror degrees of freedom. 
In the case of Eq. (|T^ we can integrate over the degree of freedom of the mirror to 



obtain the marginal distribution of the field mode as ||T0| 

2 



P{X) 



j7r/4| 



«o) 



1 

P+{X) + P_(X) + 2^P+(X)P_(X)sin (4X|ao| sin (arg 



ao) 



(23) 



.1/2 




1 exp 


-2X - |ao| 



. *2^ 



while in the case of Eq. (|T8D the marginal distribution for the field mode is in effect 
a conditional probability 

2 



P{X\yt' 



(24) 



whose explicit expression is the same as in Eq. (^31), apart an extra phase factor in 
the coherent sate which gives the interference pattern along a direction depending on 
the result of the measurement as well. 



5 Damped Mode Equation and Solutions 

Let us now consider the proposed model as an open system interacting with the 



"rest of Universe" p2[. We will study only the case in which the radiation mode 
relaxes much faster than the mirror (the opposite case, i.e. the mirror that relaxes 
much faster than the cavity mode, does not show any quantum features due to the 
thermalization effects). Moreover, since, in order to see the Schrodinger cats, we are 
interested to short time behaviour (i.e. times much shorter than the typical radiation 
relaxation time), we can consider the mirror practically not affected by any damping. 
Hence, the master equation for the whole system will be taken in the form 

p='-[p,H]+x{p), (25) 



where 23 



x(p) = ■^[2apa^ — a^ap — pa^a] , (26) 

and where we have considered the number of thermal photons to be negligible at 
optical frequencies. In our model, the damping constant 7 takes into account the 
loss of photons through the fixed mirror, so it is related to its transmissivity Tr by 
the relation 7 = cTrjIL with c the speed of light. However, since we are using a 
scaled time we should replace 7/07^ — > 7- Now, the undamped system is an exact 
solvable system with the free evolution operator t/(t) given by Eq. (|^) and obeying 
the equation lUii) = HU{t). Then, introducing a new density operator R, in a form 
similar to the interaction representation, i.e. p = URW , we may rewrite Eq. (P5|) as 



R = U^xiURW)U = x{R) , (27) 

where the operator x{R) is obtained by the following recepie: all the additional 
operators aj in the initial operators x{p) ^-^e replaced by Sj = WaiU, while the 
operator p is replaced by R. We could write down the solution of the Eq. (|27|) in the 
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form R = Rq + Y, where Rq is a constant operator, i.e. Rq = 0, and the operator Y 



satisfies the equation corresponding to ( KH) 



Y = xiY + R,) . 



{2i 



The operator _Ro represents to the free solution of the initial Eq. (p5[), i.e. without the 
term x{p)- Till now we only rewrote the master equation in another representation 
and it is still an exact equation. However, Eq. ( p7D is appropriate to apply the Born 
iteration procedure provided that the the damping term x{p) is small enough 
to be considered as a perturbative one (this could be the case since the parameter 
7 has to be small in order to achieve the Schrodinger cats). Then we could try to 
solve Eq. (|28|) simply by replacing in the r.h.s. the operator x(F + Rq) by x(-Ro); 
i.e. performing the first Born approximation. The solution is immediate, and the 
operator R assumes the form 

ft 



R{t) = Ro+ f x{Ro, r)dT . 
Jo 

It means that the initial density operator p{t) becomes 



(29) 



Pit) = Po{t) + p^{t) 



(30) 



where the term po{t) is the density operator of the free motion 

Po(t) = f/(t)po(0)f/^(t) 



(31) 



with initial density matrix po(0) = p(0). The correction term p^{t) has the form 



Pyit) = U{t) 



x{Ro,T)dr 



f/t(t) 



(32) 



or more explicitely 



a)apQ{t) + pQ{t)a)a 



(33) 



The range of validity of the above approximation is determined by the requierement 
p^(t) << po{t). Below, it will become more clear that it works for 7|aoP^ << 1- It 
is also easy to check that Tr{p^} = 0, then p{t) is always normalized to unity. Let us 
now try to find the marginal distribution at the particular times t* and t' discussed 
in Sec. ^. By means of Eqs. (|33|), (|22D and (p!4D, after lenghtly but straithforward 



algebra, one obtains 



P{X) = {X\Tr^{poin + p.imx) 

2 \ I o °° 9^(p+'J)/2 

- e-l"«l y -—HJV2X)HJV2X)e''''^''''^P-''^ 

where Hp are the Hermite polynomials. 



(34) 



A 



\P+(1 



(35) 



and finally 



Ip,iit*)= f c^re-*P^(**-")l(P-^). (36) 



In Eq. (0) the first term inside the curly brackets comes from po and is related to the 
undamped motion, while the other is the perturbative term due to the environmental 
coupling. Due to the fact that at the times t* the two subsystems (i.e. radiation 
cavity mode and mirror) are disentangled, the thermal effects do not destroy the cat 
state as can be seen in the above equations. The decoherence depends only on the 
leakage of photons through the fixed mirror. 

In the case of cat states generated by conditional measurement the expression for 
the conditional probability in presence of damping has almost the same structure of 
Eq. (|53D, and can be obtained by using Eqs. (|5B|), ( P2| ) and (|i 

P{X\y,) = {X\{y,\poit') + p,it')\y,)\X) 
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c'i- 

.IT 



2s 1/2 



°o 2-{p+g)/'2 



E 



p,q=0 



plql 



x\c^or'\^{yt'\PT\yt') + l 



Ap^gIp^g{t')\ao\'^ - Aq^p——t'{yt'\pT\yt') 



(37) 



where 



dr exp{-t[2E{t' - r) + F{t')F{t' - r) sm(r/2)](j9 - q)} 
X (yt, -F(t'-r)sin(r/2)|pTbt'-F(t'-r)sin(r/2)), (38) 



and, due to Eq. (||), the following general expression holds |p5 

•2 1 - 



(3^Ipt|3^)= (^)'(i 



71 1 + Z 



exp 



-23^' 



(39) 



C is a constant needed for the normalization after the projection in the measurement 
process, and it can be obtained by performing the integration over the X variable of 
Eq. (|37D with the aid of the completness formula for the Hermite polynomials 



C = {{yt'\PT\yt') + i\ao\^ [ip,g=pit') - {yt'\PT\yt')t']} ' . (40) 

It is easy to note that the correction term in both solutions (|3^) and (|37|) remains 
smaller than the undamped term provided 7|Q;opt << 1. Equation (|37|) shows a 
dependence of the decoherence effects also on the thermal state of the mirror (i.e. its 
temperature). In Figs. (1) and (2), we show respectively P{X) and P{X\yt' = 0) 
(solid lines) of Eqs. ( |3^ and ( p7| ) contrasted with the same in absence of damping 
(dashed lines). We may see that in the case of the cat state created at t* = 27i, i.e. Fig. 
(1), the coherence has been almost totally washed out, due to the long time needed 
for the formation; while the conditional measurement could be used to generate the 
superposition at shorter time, in Fig. (2) t' = 3tt/2, preserving the coherence effects. 
In this case, however, one should pay attention to the thermal effect of the mirror. 
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To this end, let us consider more closely the case of = 0, which is a high probable 
value for the mirror quadrature measurement. The normalization factor on the r.h.s. 
of Eq. (|39D is a common factor that can be eliminated in Eq. ( pTf ) by using Eq. 
(|40|), while the exponentioal factor remains in the integral of Eq. (|38|) only. As z 
approaches the value 1, i.e. the temperature increases, it tends to becomes unity. 
This means that the thermal effects tend to destroy the coherence only up to a value 
of temperature, above which the interference fringes become insensitive (dotted line 
of Fig. (2)). Of course analogous discussions can be made for other values of the 
mirror quadrature yt'- 

We also note from both Fig. (1) and Fig. (2) that, as the dissipation becomes 
relevant, two gaussian peaks centered around the mean number of photons, and 
which are typical of the orthogonal quadrature, appear. This is essentially due to the 
rotation in the phase space introduced by the damping term p^. In fact, as can be 
seen in Eq. (^), it involves an integration over the time which leads to a distribution 
whose contributions come from various field phases, i.e. from different quadratures. 

6 Detection of Quantum Coherence 

In this section we will show that the above discussed model could also be used to 
reveal the quantum coherence. 



According to Ref. the photon number statistics of the radiation field could 
be opportunely used as signature of the presence of Schrodinger cat states. On the 
other hand, in the presented model, a measurement of the mirror's momentum p 
allows us to get the photon number statistics in an indirect way [^. In particular 
the signal could be represented by the number a^a of photons of the radiation mode, 
and the meter by the momentum of the movable mirror; the out of phase quadrature 
coupled to the photon number (see Eq. (H)). 
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Our purpose should be to detect the Schrodinger cat immediately after its genera- 
tion inside the cavity, at time t* (or t' if one uses conditional measurement generation); 
nevertheless in both cases the two subsystems, i.e. the mirror and the radiation mode, 
are disentangled (as can be seen in Eqs. (|14D and (p!8D ), so no information can be 
extracted in indirect way. Then, we must address the measurement to get something 
which is slightly different from the Schrodinger cat state, but still having quantum 
coherence features. To this end, let us consider at first the entanglement between the 
signal and the meter, which could be described by the correlation function defined as 
follow ^ 

_ \{a^ap) - (a^a)(p)P , . 

where V means the variance. This quantity shows how good is the scheme as a 
measurement device, and should be equal to one for a perfect scheme. By performing 
the expetaction values using Eqs. (|33|) and (^ we obtain 



2 






sin^ t + 47 


sin^ (1) 


sint 




l + Nth + 2\ao 


"^K^ sin^ t 


(l-7t) + 


"0 


2K2 2[2t-8sint + 3sin(2t)] 



(42) 



Thus Cs^m is a function of t depending also on k, which is a constant that contains 
all the external parameters. Fig. (3) illustrates the typical behaviour of Cs,m versus 
t, showing the effects of dissipation as well as the thermal ones. From this figure it 
is obviously that higher values of Cs,m for times closer to 0, tt, 2ti could be achieved 
by increasing the value of k or laol, but we must take into account that the number 
of photons plays a delicate role in the dissipation effect. 

Let us now consider a time at which the radiation is entangled with the mirror, 
then its state, in absence of loss, by Eq. (^), will be 

pt'\t) = Tr^{po(t)} 
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and furthermore if E{t) satisfies the condition (|T2|) for that time, it becomes 
Po'^^^it) = \ I dyt{yt\pT\yt) 



2 

X fe"*5 



(44) 

which represents not a "pure" cat state, but that one whose phase is still convoluted 
with the mirror motion and at which we may refer as "pseudo-cat" state. This 
latter, however, has the advantage of being detected, since it does not imply any 
disentanglement. It is worth to remark that the dephasing effect due to the factor 
exp{iFyt), which degrades the pure cat into a pseudo-cat state, is considerable only 
for those values of yt contained under the gaussian state of the mirror. Then the 
temperature can emphasize this negative effect, since it introduces highest mirror 
number states, i.e. gaussians with larger width. On the other hand, in order to 
reduce this effect, it is also preferable to have the smallest possible values of F{t). 
These are accessible only at times near to 27r (see eq. (^). Thus, in order to realize 
the measurement, the choice of the mesurement time t and the value of k should be 
made to fulfill simultaneously the following requirements: Eq. (p!2|), the highest value 
of Cg^m, and the smallest value of F. Of course the detection should be performed 
at time much shorter than the typical cavity lifetime 7"^ = 2L/cTr, but also longer 
than the photon cavity fly time 2L/c, to ensure the presence of photons inside the 
cavity. 

Let us now suppose to have found the desired t and k, then we revise the mea- 
surement strategy of Ref. IPB] for the detection of quantum macroscopic coherence. 



A coherent field la,.), the "reference", is added to the pseudo-cat state, immedi- 
ately before the measurement, so that the resulting field in the cavity at the time of 
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measurement is 



1 




(45) 



Af 



where D is the displacement operator and Af is a normahzation constant. 

After the injection of the reference field, the photon number distribution in the 
cavity becomes 



where 0(7) indicates the perturbative terms proportional to the first power of 7 that 
we have omitted for space reasons. 

Let us now consider separately two cases. When ao and ar have the same phase 
the photon distribution, denoted by Vin{n), as consequence of the first term in Eq. 
(|46|), which is the dominant one, should appear as the sum of two quasi Poissonian 
distributions peaked around n = lao + c^rP and n = | — ao + «rP, with the tails due to 
the smearing effect of the gaussian integral. In fact, in Eq. (^) , the interference part 
will be negligible provide to have la^l >> 1. An interesting situation arises when 
and ar have the same amplitude, then 



In that case, neglecting the perturbation terms, Vin{n) consists of a very sharp distri- 
bution centered at n = 0, which is a (5-like peak for a pure cat state, and a distribution 
peaked around n = 4|aoP- The existence of two separate peaks in the in-phase sum 
field is the proof of the existence of two classical fields within the cavity. However, 




(46) 




where 



c± = 2 ± 2 cos {F{t)yt) . 



(48) 
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it does not prove that these two fields are in a coherent quantum mechanical super- 
position. So we need to consider also the case when ao and ctr are 7i/2 out of phase, 
for which we have 

Voutin) = i^^^— / dyt{yt\pT\yt) 



N2 



+ s:!e-l"«l'^- + 24/V/'e-2|"«l'3f?{-z(-2)"}l I + 0(7) (49) 



where now 

s± = 2±2sin(F(%i) . (50) 



In this case the interference in the term in Eq. (^6|) becomes important, in fact 
Voutin)., again neglecting the perturbation terms, exhibits a Poisson envelope with 
strong oscillations, signaling the coherence effect. The above discussed dephasing 
effect in the pseudo-cat tends to wash out the oscillations and to transform the 
Poisson envelope in a gaussian one. Of course in both cases [in and out) also the 
damping terms cause a degradation of the signal. 

In Fig. (4) we show Vin{n) for a pseudo-cat state a) which resembles that one for a 
pure cat state, contrasted with the same in presence of damping at zero temperature 
b) and at finite temperature c). Fig. (5) illustrates the same situations for Vout{n). 
Both figures are obtained using t = 0.84 x 2tt and k = 0.5 for which one has F = 0.48 
and Cs^m = 0.85 (at zero temperature, while it is reduced to 0.55 when Nth = 2). 

The Viniji) and Voutiji) distributions can actually be measured detecting the mo- 
mentum of the mirror, of course the measurement process is destructive, hence the 
state has to be reprepared for each measurement, and a large number of measure- 
ments should be performed to reach the desired statistics. Then, from these output 
distributions, one can recognize a signature of quantum coherence as in Fig. (4) 
and Fig. (5), provided to have small dissipation and very low temperature, which is 
needed also to guarantee a sufficient signal meter correlation (Fig. (3)). 
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Finally, to effectively visualize the presence of interference fringes in the phase 
space, we would consider the marginal distribution for the pseudo-cat. This proba- 
bility, obtained through the expectation value {X\TTm{poit) + P'y{'t)}\^) and using 
Eqs.(ll) and (H), will be 

/ 2 \ 1/2 oo 2-(P+'?)/2 

PP^(X)= - e-l"«l'-2^' y -— if„(V2X)iy„(y2X)e*^''S"n(P-'?) 



\p+g\^ , 7 



Ap,gIp,g{t)\ao\ -A, 



2 A P+Q, 



q,p 



exp 



F\t)ip-q)\l + z) 



(51) 



where, with the superscript pc we refer to the pseudo-cat state. It is clear from the 
last exponential factor how the thermal phonons of the mirror tend to rapidly destroy 
the coherence effect. 

In Fig. (6) we show the marginal distribution P'P'^{X) of Eq. ( |5TD for various 
situations, using the above discussed values of parameters i.e. t = 0.84 x 27r, n = 0.5. 
From this picture we may note that the interference pattern of the pseudo-cat state is 
almost the same of the pure one and is still preserved at the time of measurement, even 
in the presence of loss provided to have a very small number of thermal excitations 
in the mechanical oscillator. 

7 Conclusion 

We have proposed the use of an optomechanical model for the generation of optical 
Schrodinger cat states. We have also presented a new scheme to reveal the quantum 
macroscopic coherence, based on the new states named pseudo-cats that could be 
intended as a sort of cat states which could be recognized before their "natural birth". 
Thus the model is substantially able to produce and also to detect interference effects 
without introducing different couplings, but one should pay attention to the different 
sources of dissipation. 
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We would also point out that the studied system could be implemented for ex- 
ample idealizing the movable mirror as a piezoelectric crystal [^. The above used 
values of t, n and 7 (in the various figures), could be reached for example with the 
following set of parameters uJc ~ 10^^ s~^, oOm ~ 10^ s~^, m ^ 10^^^ Kg, L ^ 1.5 m, 
Tr ^ 10^^ and T ^ 10~^K. Of course, other choices satisfying the above mentioned 
criteria can be made giving the same qualitative results. We are aware that a delicate 
point could be the realization of the mechanical oscillator with a very small mass, but 
we would remark that the mass parameter could also be interpreted as an effective 
value coming from the density of the vibrational modes of the mechanical oscillator 



30| . Furthermore, the discussed model could be improved by inserting an active 
Kerr medium inside a cavity which enhances the nonlinear effects, slowing down the 
decoherence. 

Finally, even if we have not coupled the system under study with an external 
readout apparatus able to measure the momentum of the moving mirror, we think 
that the presented model represents an interesting alternative way to approach, also 
in the experimental sense, the quantum macroscopic coherence phenomena. 
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FIGURE CAPTIONS 

Fig. 1 The marginal distribution P{X) is plotted as function of the quadrature 
variable X ior k = 0.5, |q;o| = \/7, at t* = 27r and in two different cases: 7 = 
(dashed hne) and 7 = 2 x 10"^ (solid hne). 

Fig. 2 The marginal distribution P{X\yt/) is plotted as function of the quadrature 
variable X for yt^ = Q, k = 0.52, |q;o| = \/7, at t' = 3tt/2 and in three different cases: 
7 = 0, A^tft = (dashed line); 7 = 2 x 10"^ Nth = (solid line); 7 = 2 x 10-^, Nth > 20 
(dotted line). 

Fig. 3 The correlation coefficient Cs,m is plotted against the time for k = 0.5 in the 
case of 7 = 0, Nth = (dashed line); 7 = 10"^ Nth = (sohd line); 7 = 10"^ Nth = 2 
(dotted line). 

Fig. 4 The distribution Vin{n) vs. the photon number is plotted for a pseudo- 
cat with K, — 0.5, IcKol = V7 and t — 0.84 x 27r in the case of 7 = 0,Nth — a); 
7 = 10-2, AT,^ = b); 7 = 10-2, Nth = 2 c). 

Fig. 5 The same of Fig. 4, but for Voutip)- 

Fig. 6 The marginal distribution P*"^(X) is plotted for k — 0.5, |q;o| = -\/7, 
t = 0.84 X 27r and the values of 7 and Nth indicated in the figure. It is also compared 
with the distribution for a pure cat. 
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